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The problem of electrostatic screening of a charged line by undoped or weakly doped graphene is 
treated beyond the linear-response theory. The induced electron density is found to be approximately 
doping independent, n(x) ~ x~‘^ log^ x, at intermediate distances x from the charged line. At 
larger x, twin p-n junctions may form if the external perturbation is repulsive for graphene charge 
carriers. The effect of such inhomogeneities on conductance and quantum capacitance of graphene 
is calculated. The results are relevant for transport properties of graphene grain boundaries and for 
local electrostatic control of graphene with ultrathin gates. 


I. INTRODUCTION 


One of the properties that make graphene an attractive 
platform for electronic devices is tunability of its charge 
carrier density through electrostatic gating. The gates 
can be brought into immediate proximity of graphene, 
which enables one to control doping of this material 
very efficiently. Local gating on ultrasmall lengthscales 
is attainable by utilizing nanowire^ or nanotube^ gates 
[Fig. 1(a)]. Physical insight into fundamental charac¬ 
teristics of such devices can be gained from a simplified 
problem of how graphene responds to a linelike external 
charge [Fig. 1(c)]. This problem is also relevant for under¬ 
standing properties of grain-boundary defects [Fig. 1(b)] 
in graphene grown by chemical-vapor deposition. 

The problem of screening of a linelike charge by elec¬ 
trons in graphene is an interesting challenge because the 
usual linear-response theory^ fails when the line is highly 
charged and/or when graphene is lightly doped. Previ¬ 
ous studies of one-dimensional charge perturbations in 
graphene eschewed solving this difficult problem resort¬ 
ing instead to ad hoc approximations for the induced 
carrier density profile.Accurate determination of this 
profile requires numerical calculations, e.g., finding the 
self-consistent solution of the Dirac equation for electron 
wavefunctions and the Poisson equation for the electro¬ 
static potential ff>(x) as a function of the in-plane coor¬ 
dinate X transverse to the charged line. However, if ff>(x) 
varies smoothly on the scale of the local Fermi length 
kp^{x), a simpler approach based on the Thomas-Fermi 
approximation (TFA) can be applied.^ We show that ap¬ 
proximate solutions of the TFA equations for the electron 
density can be derived in certain limits. Using these ana¬ 
lytical solutions and numerical simulations, we compute 
two other important observables amenable to experimen¬ 
tal probes: the conductance and the gate capacitance. 

Let us introduce our key notations and assumptions. 
We denote the linear charge density of the external per¬ 
turbation by eX. Without loss of generality, we take A to 
be positive (unless specified otherwise). We assume that 
the unperturbed electron density Uoo of graphene is uni¬ 
form. To distinguish between n-type and p-type doping, 
we define the Fermi momentum corresponding to rioo as 


a signed quantity, 

koo = sgn(noo)|7rnoo|^/^. (1) 

We assume that the external charge enext(^) has a 
Lorentzian density distribution, 

/ \ X a . . 

^ext(^) = -^--2- 

TT x^ P 

The actual profile of the external charge may of course be 
somewhat different. For example, the charge distribution 
of a grain boundary probably does not have power-law 
tails. However, the role of parameter a in Eq. (2) is 


(a) nanowire gate 



FIG. 1. (Color online) Models of linear charged perturba¬ 
tions in graphene devices, (a) Nanowire top gate. The bot¬ 
tom gate separated from graphene by an insulator of high 
dielectric constant k may be useful for additional control, (b) 
Grain boundary (pentagon-heptagon chain) with charged ad¬ 
sorbates (circles), (c) Charged string off the graphene plane. 
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mainly to regularize the response of graphene at very 
short distances,^ i.e., it serves as a short-distance cutoff. 
As long as we are not interested in microscopic physics 
at \x\ < a, Eq. (2) can be adopted as a convenient model. 
In all examples mentioned above (nanowire and nanotube 
gates and also grain boundaries in graphene) realistically 
achievable a can be as small as a few nanometers. The 
particular functional form of Eq. (2) corresponds to the 
idealized model shown in Eig. 1(c) where the external 
charge is located off the graphene plane and is truly one¬ 
dimensional. This can be seen from the fact that the 
electrostatic potential created in the graphene plane by 
the out-of-plane charged line is equal to that created by 
the in-plane Lorentzian charge distribution (2): 




e f 1 

tix) = - j dx'n^^tix')\og (3) 

(4) 


e\ ^ 1 

= ^ —2 • 


In general, the effective width parameter a should be 
taken as the larger of the actual physical width of the 
charge distribution and its distance to the graphene 
plane. In this article we assume that electron-electron 
interaction in graphene is weak, i.e., we consider the di¬ 
mensionless coupling constant a = jnhv a small pa¬ 
rameter in the problem.Here v = lO^crns”^ is the 
graphene Eermi velocity. By choosing a substrate with 
a large dielectric constant n [Eig. 1(a)], it is possible to 
make a 2.2/n significantly smaller than unity. How¬ 
ever, it is difficult to make a truly small, so as a rule 
we do not treat log a as a small parameter. Einally, we 
assume that the graphene is not too highly doped. 


l^ool « 


1 

aa 


(5) 


Depending on the relation between A, and x, the 
response of graphene can be either weak or strong and 
either linear or nonlinear (Eig. 2). The degree of nonlin¬ 
earity is controlled by the dimensionless parameter 


~ a\ 

^ ~ k ’ 


(6) 


Linear screening is realized if A ^ 1 (the bottom part of 
the ‘weak’ region in Eig. 2) where the induced electron 
density 


^ind(^) = ri{x) - 


(7) 


scales linearly with A. On the other hand, if A > 1, 
the response is nonlinear. A conspicuous manifesta¬ 
tion of the nonlinearity is found in the region labeled 
‘strong’ in Eig. 2, where the induced density can be ap¬ 
proximated by a ‘universal’ (doping-independent) form 
^ind(^) ^ x~^log^ X. In the large-x ‘asymptotic’ regime 
of Eig. 2, the induced density exhibit a power-law decay 



FIG. 2. Schematic diagram of screening regimes. In the 
‘weak’ region screening of the external potential is poor; the 
response is linear [Eq. (15)] at A < koo/a and nonlinear 
[Eq. (24)] at larger A. In the ‘strong’ regime the external 
potential is greatly reduced and the induced density profile 
is given by Eq. (24). In the ‘asymptotic’ regime the den¬ 
sity profile follows Eq. (8). In the ‘perfect’ regime graphene 
maintains local charge neutrality, apart from small corrections 
[Eq. (30)]. This diagram is drawn assuming graphene is not 
too heavily doped, Eq. (5); otherwise, the ‘weak’ and ‘strong’ 
screening regions would disappear. 


where the dependence of 6 on A is linear [Eq. (11) of 
Sec. H] if A is small and logarithmic if A is large [Eq. (27)]. 
The ‘perfect’ screening region is included in Eig. 2 for 
completeness. Here graphene maintains charge neutrality 
locally, i.e., the induced density is close to the external 
one, nind(^) — ^ext (^)5 so that the model assumption (2) 
must be critically revisited. 

The crossovers among the predicted screening regimes 
can be systematically studied in experiments using de¬ 
vices that have both top and bottom gates, Eig. 1(a). 
High linear charge densities with aX ^ (lnm)“^ are 
quite feasible to achieve with top gate voltages E ^ 1V. 
Eor lightly doped graphene, n^o = 10^^ cm“^, the corre¬ 
sponding A ^ 20 is deep in the nonlinear regime. 

In the remainder of this article we derive detailed for¬ 
mulas for the carrier density profiles, verify them by nu¬ 
merical simulations, then make predictions for capaci¬ 
tance and transport measurements. 


II. LINEAR RESPONSE 


Linear screening of linelike charges by doped graphene 
has been studied in previous literature.^ We include a 
brief summary of the relevant results for later comparison 
with our nonlinear response theory. Linear screening is 
realized when the external charge is not too high or when 
graphene doping level is not too low. The quantitative 
criterion A ^ 1 is derived below. Within the linear- 
response theory, the induced electron density is given by 


oo 

'^ind(^) ~ J 


dq 

27r 


1 - 
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FIG. 3. (Color online) Comparison of induced electron density profile from the numerical solutions of the TFA equations (blue 
dots) and analytical formulas (black lines) in three screening regimes. The red curves (labeled ‘Perfect’) in each panel represent 
ne^t{x). Units for density n, distance x, and linear number density A are respectively (dvr^a^a^)”^, a, and (47r^a^a)~^. (a) 
Linear screening regime realized for A = 0.1 and Uoo = 4. The induced density is much smaller than rioo- (b) Nonlinear regime 
realized for A = 0.4 and rioo = 0. Density profile characterized by screening length Xg ~ 25. (c) Near-perfect screening for 
A = 1000 and rioo = 0. Note the double logarithmic scale. 


where the term is the Fourier transform of the 

effective external charge next(^) [Eq. (2)] at momentum 
q. As we are primarily concerned with distances x ^ 
where the TFA is valid, the dielectric function of 
graphene can be approximated by^^“^^ 


Refinement of Eq. (10) can be obtained through the 
random-phase approximation (RPA). Within the RPA, 
the dielectric function of graphene coincides with Eq. (10) 
at 0 < jg'l < 2 \ koo\ but at \q\ > 2 \ koo\ it is given by a 
different formula^ 


e(g) = l + ^, qxF = ^cx\koo\‘ (10) 

A particularly simple analytical expression for nind(^) 
can be derived in the ‘asymptotic’ regime, cf. Eig. 2. 
Carrying out the integration in Eq. (9) by the steepest- 
descent method, one finds the leading-order approxima¬ 
tion for the induced density to be^ 

nind(a:)~^, b=^{a + q:^p). (11) 

(In contrast, nind oc x~^ for a pointlike charge perturba¬ 
tion.^^) Note that the coefficient b in Eq. (11) is much 
larger than Xa/n under the assumed condition (5), so 
that nind(^) ^ ^ext(^)- Such an ‘overshoot’ is typical for 
screening of localized perturbations in metals. Metallic 
systems possess the overall charge neutrality. However, 
at distances shorter than the local screening length from 
the perturbation screening is necessarily weak. There¬ 
fore, there is a missing charge at short distances, which 
has to be compensated at large x. However, the electric 
field is not overscreened: it is of the same sign as the 
external one but reduced in magnitude. 

The analytical Eq. (11) agrees well with our numerical 
simulations shown in Eig. 3(a). Eor these simulations we 
used previously developed codes^^ with suitable modifi¬ 
cations. In brief, the electron density n{x) to be found 
was defined on a grid of x with periodic boundary condi¬ 
tions. The solution was obtained by minimizing the total 
energy of the system (kinetic plus electrostatic) within 
the TEA using standard technical computing software. 


/ X _ 1 , ^TF 



+ a cos 


2feoo 

Q 


\q\ > ‘^\koc \. 


( 12 ) 


Notably, e 1 + 7ra/2 becomes doping independent at 
|g| ^ 2 |/coo| where the response is dominated by virtual 
interband transitions. (Eor corrections to the last result 
beyond RPA, see Refs. 10 and 15.) 

Substituting Eq. (12) in Eq. (9) and using contour in¬ 
tegration techniques to evaluate the integral, we find the 
RPA correction to nind- 


A RPA / \ 7 cos (2 I ffi 4) 17 I 

- h --[^75- - , \kocx\ » 1: 

\kooX\^^ 


bi = 


A a\ko 


20f (l + 2a)2 ’ 


which is a particular case of the Eriedel oscillations, 
intermediate distances, Eq. (9) yields^ 


(13) 

(14) 
16 At 


nind (a:) 


A(?tf 

TT 


log 


0.561 

IgxF^I 


1 < lA^oo^l <C (a ^. (15) 


Einally, let us estimate the region of validity of the linear- 
response theory. This theory applies if the induced car¬ 
rier density is smaller than the original one, nind ^ ^00 
or, equivalently, if the local Eermi momentum, 

kpix) = sgn(n(x)) |7rn(x)|^/^ 


(16) 
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is perturbed slightly compared to koo [Eq. (1)]. Naively, 
one may require the condition [kpix) — koo]/koo ^ 1 to 
hold at ah distances of interest, x ^ a. In fact, the 
validity region is wider because the nonlinearity affects 
only the response at momenta q <2 |A:oo|- (As mentioned 
above, the response at q > 2\koo\ is essentially doping- 
independent.) Therefore, the smallest number we should 
use in the argument of nind(^) for out estimate is x ^ 
|/Coo|~^- From Eqs. (7), (15), and (16) we get 


max 


krjx) _ ^ 
koo 


2aX 

koo 


log 



(17) 


Neglecting the logarithmic factor, which is never large in 
practice, we arrive at |A| ^ 1 as the criterion of linear 
screening. 


III. NONLINEAR RESPONSE 


In this Section we treat a more difficult case A ^ 1 
where screening is nonlinear. Our approach to this prob¬ 
lem is to solve the TEA equations analytically and nu¬ 
merically. The hrst of these equations is 

/i(x) — e4>(x) = 0 , (18) 

where 


/i(x) = hvkpix) (19) 

is the local chemical potential of graphene (assuming the 
linear Dirac dispersion of quasiparticles). The second 
equation links the total charge density and the electro¬ 
static potential, 

2e f 

$(x) = — dx'log \x - a:'|[ni„d(a;') - ne^t{x')]. (20) 

This relation can be inverted by exploiting techniques 
from the theory of analytic functions: 

CX3 

/ \ \ ^ ^ f x'dx' . 

'^ind(^) ~ '^ext(^) “ ^ ^ 2 ’ (^^) 

TT^e J x’^ — x^ dx' 

0 

where V stands for the Cauchy principal value. If desired, 
Eqs. (18)-(21) can be combined into a single nonlinear 
integral equation for n{x). The TEA is valid if ^ 

^k-p\x)^l. (22) 

dx 

The problem we want to solve can be separated into two 
parts, depending on whether nind(^) is greater or smaller 
than rioo- The latter situation occurs at large x, where 
we expect screening behavior akin to linear response. In¬ 
deed, it is easily seen that nind(^) follows Eq. (8) provided 
the integral in Eq. (21) is dominated by x' x. Invok¬ 
ing Eq. (18), we then obtain the asymptotic behavior 
4> ^ x~‘^ for X ^ Xoo, where 

Xoo = (23) 


and b is to be determined below. 

The analytical form of nind{x) at x <C Xoo is not im¬ 
mediately obvious. Fortunately, were are able to hnd (by 
trial and error) the following approximate solution: 


^ind (^) 


1 1 

47r^(a^ {x — Xg)^ 


log^ 



(24) 


which is characterized by a nonlinear screening length Xg . 
This length is found from the argument that at small x 
the external held is nearly unscreened, so that the to¬ 
tal and the external potentials differ only by some con¬ 
stant: 4>(x) 4>ext((^) + const = —2X{e‘^/K) logx + const 

[Eq. (4)]. Comparing with Eq. (24) at x <C we get 

1 

^ 47ra‘^X 

From Eqs. (22) and (25) we see that Eq. (24) is valid 
ut X 271 ax g. At smaller x the density and the Fermi 
momentum presumably tend to a hnite maximum, i.e. 


lUdixkFix) 


kpi^^TraXg) 7^ 2(aAlog 



(26) 


The last logarithmic factor is valid if a 1; otherwise, 
it should be replaced by a number of the order of unity. 
Equation (26) is consistent with Eq. (17) for the linear 
regime. Hence, in both regimes the nonlinearity parame¬ 
ter A [Eq. (6)] has the physical meaning of the maximum 
relative change in the Fermi momentum {kp — koo)/koo 
caused by the perturbation. 

To hx the so far undetermined coefficient 5, we require 
a smooth matching of Eq. (8) and Eq. (24), which yields 

lib1^1' = 

More accurate expression for n[nd{x) can be obtained 
by iterations using Eq. (24) as the input to the TEA equa¬ 
tions. Namely, substituting it into Eq. (18), we can get 
4>(x), which we can then insert into Eq. (21) to obtain an 
improved approximation for n[nd{x). The hrst iteration 
yields 




(x) 


1 

4:7r^a‘^ 


TT^ 1 
2 (x-hXs)^ 


+ log^ 



b +a;^ _ logObab 

(x^ — x^)^ x‘^ — x‘1 


(28) 


which demonstrates a good agreement with our numeri¬ 
cal simulations, see Fig. 3(b). 

For completeness, we consider the case of 


dTTQf^a ’ 

where nonlinear screening affects distances shorter than 
our cutoff length a. At such large A, a highly doped re¬ 
gion appears at x < a, where screening is nearly perfect. 
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i.e., where the difference between nind(^) and next(^) is 
relatively small and the system is able to maintain the lo¬ 
cal charge neutrality. The solution for nind(^) can again 
be obtained by iterations. The input to Eq. (21) is now 
^ind(^) = ^ext(^)- T^c first corrcctiou is 


^ind(^) ^ext(^) 


^/Xa X sinh ^ ^ 
(x^ + a^)i 


1 

x‘^ + 


. (30) 


Its range of validity x < Xp can be estimated from the 
requirement that ai x = Xp this correction is no longer 
small, which implies 


\og{xp/a) 


^ ira 


\/Aa. 


(31) 


Furthermore, one can check that at x = Xp Eq. (30) 
matches by the order of magnitude with nind(^) given by 
Eq. (24) with Xs set to a. This leads to us conclude that 
at X > Xp Eq. (24) and then at x > x^o Eq- (8) must 
still hold with Xg ^ a. Because of the exponentially large 
magnitude of we could verify numerically only the 
X < Xp case. Comparison of simulations with Eq. (30) in 
Fig. 3(c) indeed shows a good agreement. 

If Xp > y/Xa/koo^ which corresponds to 


A > 


1 

TT^Qf^a 


log" 



(32) 


the screening should be nearly ‘perfect’ for all x > a, 
cf. Fig. 2. 

When A and rioo have opposite signs, twinp-n junctions 
appear at \x\ ^ Xoo- The density profile near a graphene 
p-n junction has been studied in Ref. 13. These results 
should more or less carry over to the present problem, so 
they will not be repeated. In the following Sections we 
focus on computing the effect of nonlinear screening on 
transport and capacitance characteristics of the system. 


IV. QUANTUM CAPACITANCE 

The charge density A is a natural parameter for 
graphene grain boundaries [Fig. 1(b)]. In contrast, 
linelike perturbation created by means of narrow gates 
[Fig. 1(a)] are controlled by the voltage V with respect 
to graphene while A has to be found by integrating the 
differential gate capacitance 

0(V) = e§. (33) 

It is important that V is not simply equal to the electro¬ 
static potential difference A4> between the gate and the 
graphene sample. It has another contribution from the 
graphene chemical potential: 

U = A4> + (p/e). (34) 

As a result, the differential gate capacitance C has two 
components, the classical (or geometrical) one Cg and 
the so-called quantum one Cg, which add in series: 

C-^=Cp + Cp. (35) 


Our goal in this Section is to derive the quantum capaci¬ 
tance of a device with an ultranarrow top gate. Fig. 1(a). 
However, it is useful to review the conventional planar- 
gate geometry first. Here the classical capacitance per 
unit area Cg = /^/(dTrd) is inversely proportional to the 
separation d between graphene and the planar gate. In 
turn, the quantum capacitance is proportional to the 
thermodynamical density of states, Cg = e‘^{dn/dij). 
This quantity can also be written in terms of the inverse 
Thomas-Fermi screening length gxF- 

Cg = ^ Qtf • (36) 

The net result of having a finite screening length is equiv¬ 
alent to replacing the physical gate-graphene separation 
by an effective one: 


deff — d + ^ ^XF ‘ 


(37) 


Although in this article we use the free-fermion approxi¬ 
mation [Eq. (10)] for the linear-response screening length 
in reality it is modified by many-body interactions 
and disorder (see, e.g.. Refs. 17 and 18). Recently quan¬ 
tum capacitance measurements have been used to probe 
such effects of graphene. 

If the gate is now a long metallic string or radius / ^ a, 
the capacitances per unit length are relevant. In order 
to derive C~^ we start with a general expression for the 
electrostatic potential difference between the string and 
the graphene sample. 


A4>(x) = ^ log 

K 




e 

H— 
n 


CXJ 

dx Tl[Yid{x ) 


log 


-h 

{x — x'Y ’ 


(38) 


which follows from Eqs. (3) and (20). If we set riind = 
next, we obtain A4> = (eA/zi:) log(4a^//^) = const. 
Hence, the geometric capacitance of the system is 


A<I> 




(39) 


which can be alternatively derived by the method of im¬ 
ages. Next, combining Eqs. (18), (20), (33), (34), (38), 
and subtracting we find 


oo 

Cp = ^ y dxn^^t{x)— ^{x). 


(40) 


We can now use this expression for analytical and nu¬ 
merical calculations. Our analytical formulas for positive 
A = X/{akoo) are as follows: 


2 


j , A 1, (41a) 

— \og{4:7ra^aX ), A 1, (41b) 
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where Ei(z) is the exponential integral. These equations 
describe, respectively, the linear and the ‘strong’ regimes 
of Fig. 2. They can be commonly written as 

(«) 

where Xsc is equal to (in the same order) and Xs . We 
conclude that the total capacitance can be modeled after 
the geometric one, 

with the effective gate-graphene separation 

aeff = a+^, (44) 

where a on the right-hand side was added by hand to 
recover the result Ueff = a expected for a perfect metal, 
Xsc = 0. The similarity of Eqs. (37) and (44) illustrates 
once again that quantum capacitance is a measure of the 
screening length of a system. Whereas in the planar-gate 
geometry this length is formally divergent for undoped 
graphene, for the linelike gate the divergence is regular¬ 
ized by nonlinearity. Numerically, we find that C~^ as a 
function of A approaches a universal envelope curve (41b) 
shown by the dashed line in Fig. 4. 

For completeness, we consider the ‘perfect’ screening 
regime where 

e$(a;) ~ next(a;) ^ (45) 

for all relevant x, so that Eq. (40) becomes 


oo 



— oo 


The inverse thermodynamic density of states dfi/dn in 
Eqs. (45) and (46) is to be evaluated at n = 'noo+^ext(^)- 
For rioo = 0 where the integrand scales as we find 

the analytical result 



27r(a\/Aa 


A> 


(47) 


which agrees with our numerical simulations (the dashed- 
dotted curve in Fig. 4). 

If A is negative, i.e., if A and k^o have opposite signs, 
the twin p-n junctions form at some A which can be esti¬ 
mated from Eq. (17) setting kpiO) to zero. This event — 
the onset of the ambipolar regime — is marked by a max¬ 
imum in which is absent in the unipolar trace for 

the same |/Coo|, see Fig 4. From dimensional arguments, 
as koo approaches zero, the height of the maximum in 
measured with respect to its plateau at A = 0 
should approach a universal number. Figure 4 suggests 
that number is about 1.5. However, such dimensional ar¬ 
guments assume the TFA is valid, which, similar to the 
case of a single p-n junction,is the case at small a. To 
treat a more typical case a 1 one needs to go beyond 
the TFA, which may be a problem for future research. 



FIG. 4. (Color online) Inverse quantum capacitance as a func¬ 
tion of the gate charge. The curves are from numerical cal¬ 
culations for koo specified in the legend. The units of A and 
koo are (47r^o^a)“^ and {7T^aa)~^. The triangles, the dashed 
line, and the dash-dotted line correspond, respectively, to the 
asymptotic limits of the linear, nonlinear, and near-perfect 
screening, Eqs. (41a), (41b), and (47). The difference be¬ 
tween the curves and the triangles at smaller koo are due to 
finite-size effects in the simulation. The peaks of the koo < 0 
curves signal the formation of the twin p-n junctions. 


V. CONDUCTANCE 

Charged linelike defects are known to significantly in¬ 
fluence electron transport in graphene. Grain bound¬ 
aries strongly reduce the sheet conductivity of large-area 
graphene,while bipolar junctions created by nanowire 
gates cause conductance oscillations.^ In this section we 
find expressions for the graphene conductance G relevant 
for both situations. 

We consider a scattering problem for a massless Dirac 
particle with initial momentum k = |/coo|(cos^, sin^) 
subject to the potential perturbation 

- e4>(x) = /i(noo) - l^{n) = hv[koo - kpix)] . (48) 

The intermediate equation follows from the TFA, 
Eq. (18). For a weak potential, we can apply the stan¬ 
dard perturbation theory to the massless Dirac equation 
to obtain the reflection coefficient 

oo 

r{0)=itdin0 j e^^^^^kF{x)dx ^ (49) 

— oo 

which is similar to the first Born approximation formula 
for the Schrodinger equation.^ The region of validity of 
this formula can be extended beyond the perturbative 
regime if in the argument of the exponential we replace 
kx by A^f(O) cos 0, the local momentum at the x = 0 point 
where the scattering potential is the ‘most’ nonanalytic.^ 
However, this is permissible only if kpix) is real at all x, 
i.e., if all points on the quasiparticle path are classically 
allowed (no quantum tunneling occurs). 

The conductance G is found by summing the transmit- 
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tances T{9) of all the ky = kao sin0 channels: 

G=^^T(0), T{e) = l-\rm". (50) 

ky 


In the absence of scattering, A = 0, the conductance is 


_ 4e" \k^\W 

h TT 


(51) 


where W is the width of the graphene sheet. Assuming 
W ^ we can compute G at A 7 ^ 0 by replacing the 
summation with the integration over 0. However, \r{0)\ 
diverges as | 6 >| tends to 7 r/ 2 , so that the first Born approx¬ 
imation cannot be used. In fact, the absolute value of the 
exact reflection coefficient must be approaching unity in¬ 
stead of diverging. We account for this by cutting off the 
integration limits at 0 where \r{0)\ ^ 1 : 


G = Go 


u 

I 


iOde. 


(52) 


is the electric field at the junction/^ with n' being the 
density gradient. Combining these equation, we obtain 


Tpn{()) = exp sin^ > ^2 ~ log^^^ |A|. (56) 

Multiple reflections of the quasiparticles in the region 
between the p-n junctions lead to the conductance oscil¬ 
lations and resonances. The net transmittance is given 
approximately by the Fabry-Perot-like formula^ 


T{0) 


TpnjO) 

1-[1-Tp„(0)]e*^ 


(57) 


where (j) is the phase acquired by a quasiparticle after one 
roundtrip between the junctions: 



(58) 


For the linear screening regime, A <C 1 [Eq. ( 6 )], we ob¬ 
tain 


Go -G 
Go 


^ log' |A| 




l|A-i + (5-i| ’ 


|A| ^ , 

q;2 < |A| < (5. 


(53) 


Accordingly, the conductance G as a function of A ex¬ 
hibits an asymmetric maximum at A = 0 drawn schemat¬ 
ically in Fig. 5. The asymmetry becomes pronounced 
when A approaches S = 1/(2 loga”^), where the screen¬ 
ing begins to cross over into the nonlinear regime. If A is 
positive, kpix) becomes essentially independent of k^^. 
Using Eqs. (24) and (26) and changing the integration 
variable x ^ Xx in Eq. (49) one can show that r{0) does 
not depend on A any more in this regime. Actually, this 
is clear from dimensional argument: since r{0) is dimen¬ 
sionless, it may not depend on A, which has the units of 
inverse length. This implies that G ceases to decrease 
with A, leveling at a plateau, (Go — G)/Go — ttS. If A is 
negative, the situation is quite different. The scattering 
potential is repulsive. At large |A| it causes the twin p-n 
junctions to appear at \x\ ^ x^o [Eq. (27)], which act 
as tunneling barriers. Here even the modified Born ap¬ 
proximation fails completely. The transmittance Tpn{0) 
of each p-n junction is given instead by^^ 

Tpn{d) = exp sin^ 0^ , (54) 

where 

(55) 


The last estimate is obtained using Eqs. (24) and (27). 
Conductance minima arise at 0 = ( 2 m + l) 7 r, where m is 
an integer. They have the magnitude Gmin /Go r\j \f 0 Ljh 2 
because the transmittance of e ach ju nction is appreciable 
only at small angles^^ d < cf. Eq. (56). Con¬ 

ductance maxima are found at = 2 m 7 r, which cor¬ 
respond to jA^^I ^ exp( 7 rv^ 2 mn). The widths of these 
maxima and the distance from one to the next increase 
exponentially as a function of m. The heights of these 
maxima approach Go, which is a manifestation of the res¬ 
onant tunneling phenomenon. In practice, observation of 
the resonant tunneling requires samples with the mean 
free path longer than the roundtrip distance 4xoo; oth¬ 
erwise, the conductance is influenced by diffusive trans¬ 
port between and across the p-n junctions.In previous 
experiments with nanowire-gated graphene devices,^ the 
conductance maxima were found to be significantly lower 
than Go and decreasing with the top gate voltage, i.e., 
|A|, presumably due to disorder scattering. 

If the conductance can also be measured in the di¬ 
rection parallel to the gate, we expect it to show a de¬ 
pendence consisting of a smooth increase with |A| with 
superimposed small oscillations due to quantization of 
the quasi-bound resonant modes. This oscillating part 
would have the same period as the Eabry-Perot oscil¬ 
lations in the transverse conductance discussed above. 
At small gate voltages, both longitudinal and transverse 
conductances are expected to show additional fine fea¬ 
tures related to the analog of the Goos-Hanchen effect in 
graphene,which is a lateral displacement of a quasipar¬ 
ticle trajectory along the p-n interface during reflection. 
This effect can be included using a more accurate equa¬ 
tion for 0 that incorporates both the path length contri¬ 
bution expressed by Eq. (58) and the phase shift of the 
reflections at the p-n interfaces. 


F = 2.5 hva^^^\n\‘^^^ 













G/Gq 


1 



FIG. 5. Reduced conductance GIGq as a function of A 
(schematically). A monotonic decrease is expected in the 
unipolar case, A > 0, while Fabry-Perot oscillations with a 
period determined by Eq. (58) should dominate for negative 
and large A where the p-n junctions form. Parameter 5 is 
defined in the text. 


VI. DISCUSSION 

In this work we considered linelike charged perturba¬ 
tions on graphene and derived analytical expressions for 
the induced density profiles in both linear and nonlin¬ 
ear screening regimes. These results were applied to the 
analysis of two types of electronic properties. The first 
one is the quantum correction to the classical capacitance 
between the narrow gate and graphene as a function of 
the top gate voltage. Measuring this quantity will be a 
direct way for observing the crossovers among different 
screening regimes and testing our predictions. For exam¬ 
ple, we showed that the divergence of the inverse quan¬ 
tum capacitance of undoped graphene predicted from the 
naive linear-response theory will be curbed by nonlinear 
screening effects. If the gate creates a strong repulsive 
potential for charge carriers in graphene, it can induce 
twin p-n junctions. Our calculations indicate that the 
onset of this ambipolar regime is signaled by a peak in 
the inverse quantum capacitance. 

The second quantity we studied is the transverse elec¬ 
trical conductance of the system. Our predictions for 
the ambipolar regime, where the conductance oscillates 
as a function of gate voltage, include formulas for the 
maxima, minima, and the oscillation period. Our theory 
holds in the ballistic transport regime, which was diffi¬ 
cult to probe in earlier experiments on such systems.^ 
We hope that modern higher-quality devices that utilize 
graphene encapsulated in boron nitride^^ and single-wall 
nanotube gates of smallest possible diameter, would en¬ 
able a systematic investigation of nonlinear screening and 
resonant tunneling phenomena we discussed. 

Conductance of graphene with charged grain bound¬ 
aries was previously studied analytically and numerically 
in Refs. 6, 26, and 27. In Ref. 6 transport properties 


were computed modeling charged grain boundaries as 
short line segments of length W ^ k^. However, the 
screening was treated assuming that the boundaries are 
infinitely long, which seems to require the opposite in¬ 
equality. These incompatible assumptions make a direct 
comparison between our analytical results for the conduc¬ 
tance difficult. As for screening, only the linear regime 
was considered in Ref. 6, and for this our results agree. 

Our findings have further implication for experiments 
using novel scanned-probe techniques. Linear charged 
defects in the form of grain boundaries have been shown 
to reflect surface plasmon polaritons^’^^’^^ and induce 
photocurrent,both of which can be imaged with 
nanoscale resolution using scanning near-held optical mi¬ 
croscopy. Scanning tunneling microscopy is another av¬ 
enue of approach to measure local density of states. 

We will apply our theory to interpretation of such mea¬ 
surements in a future work.^ 

Finally, let us comment on the proposals^^ that 
graphene is a condensed-matter laboratory for exotic ef¬ 
fects predicted in other helds of physics. For example, 
electronic response of graphene to a pointlike charge has 
an interesting analogy to the atomic collapse of super¬ 
heavy elements.^^’^^ It has been shown^’^^’^^“^^ that in 
graphene subcritical Z < and supercritical Z > Z^ 
charges produce qualitatively different behavior of the 
screened electrostatic potential at large distances from 
the perturbation, the critical charge Zc being of the or¬ 
der oil/a. Characteristic oscillations of the local density 
of states that appear in the supercritical case have been 
recently detected experimentally.^^ Phenomena similar 
to atomic collapse have been also studied theoretically 
in the context of narrow-band gap semiconductors and 
Weyl semimetals.In turn, our problem of screening of 
a linelike charge perturbation in graphene have interest¬ 
ing analogies in cosmology (screening of a hypothetical 
cosmic string by vacuum polarization^^) and polyelec¬ 
trolyte physics (Onsager-Manning condensation of coun- 
terions^^’^^). Our results imply that in graphene non¬ 
linear screening plays a greater role for linelike charges 
compared to the pointlike ones: the former are always 
supercritical, e.g., there is no threshold A for the appear¬ 
ance of Friedel oscillations. Finally, our analytical formu¬ 
las assume a ^ 1, which can be realized using a high-/^ 
dielectric substrate [Fig. 1(a)], such as SrTiOs [Refs. 46- 
49] . For such gate dielectrics it may be important to con¬ 
sider electric-field dependence of k in nonlinear screening 
regimes. 
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